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ABSTRACT 


Ik  computer  program  was  developed  which,  is  particularly  well  suited  for 
use  k ^n ^parametric  flutter  investigations  where  a  large  quantity  of  numerical 
data  sSt*’ required.  The  underlying  analysis  used  exact,  linearized,  three- 
dimensional  aerodynamic  theory  so  the  program  can  he  used  to  investigate  the 
critical  low  supersonic  regime.  Finite  panel  arrays  can  he  treated,  and  wind 
tunnel  installations  can  he  simulated.  A  limited  quantity  of  informative 
numerical  data  is  presented.  Most  of  the  data  pertains  to  the  flutter  of  sin¬ 
gle  panels  with  free  side  edges,  with  and  without  adjacent  vertical  side  walls 
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=  modulus  of  elasticity  of  panel  material 
=  matrix  defined  in  Eq.  (38) 


v  - 


F(u)  =  function  defined  by  Eq.  (29) 

-  deflection  functions  (Eqs . (19,20) ) 
m  ^  =  integral  defined  by  Eq.  (36) 

g  =  structural  damping  coefficient 

Hg  m  ^  =  integral  defined  by  Eq.  (36) 

h  =  panel  thickness  (dimensional) 

Ijj  mu  =  aerodynamic  integral  (Eq,  (35)) 

Jm,m  =  structural  integral  (Eq.  (9)) 

K  =  kM/(32 

m  =  structural  integral  (Eq.  (9)) 


k 

L 


£,I 

M 


m,m 


N 

Pu 

P 

P 

^u;m 


Q 


q. 
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=  torque  per  unit  length  exerted  by  stringer  on  panel  edge 
=  torque  per  unit  length  exerted  by  rib  on  panel  edge 
=  time 

=  free  stream  velocity 
=  index  in  Fourier  expansion  (Eq.  (23)) 
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v  =  Poisson's  ratio 

\ 

§  =  nondimens ional  chordwise  coordinate 

p  =  mass  density  in  undisturbed  stream 
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as  chordwise  deflection  function 

0  =  displaced  coordinate  in  Fourier  expansion  of  spanwise  deflection 

shape 

as  spanwise  deflection  function 
cp  =  local  coordinate 

03  =  dimensional  flutter  frequency 


-  viii 


I.  INTRODUCTION 


There  is  at  the  present  time  a  serious  need  for  theoretically  derived 
quantitative  data  on  the  flutter  of  flat  panels  in  a  supersonic  air  stream. 
Aircraft  and  missile  designers  need  the  information  as  do  persons  interested 
in  refining  the  theory  by  comparison  with  experimental  results.  To  be  gen¬ 
erally  useful  the  data  must  apply  to  structural  configurations  of  practical 
interest,  and  must  be  extensive  enough  to  differentiate  the  effects  of  the 
numerous  nondimensional  parameters  which  characterize  the  flutter  phenomenon. 

The  flutter  problem  has  been  adequately  formulated  for  some  time  but 
limited  numerical  data  are  available  because  of  the  considerable  effort  re¬ 
quired  to  obtain  solutions  to  the  flutter  equations.  In  the  present  study  a 
computer  program  was  developed  which  obtains  solutions  to  the  flutter  equations 
with  relative  ease.  The  program  is  written  for  the  IBM  7094  computer,  and  has 
two  significant  features: 

1.  The  underlying  analysis  is  broad  in  scope  such  that  a  wide 
variety  of  physical  conditions  is  covered.  (For  example,  typical  experimental 
wind  tunnel  situations  can  be  analyzed. ) 

2.  A  technique  for  the  solution  of  the  flutter  equations  is  used 
which  is  particularly  well  suited  for  parametric  studies,  i.e.,  extensive 
parameter  variations  can  be  accomplished  with  minimum  computer  effort. 

These  two  features  make  the  computer  program  of  considerable  practical  value 
in  view  of  the  preceding  discussion. 

Computations  were  carried  out  for  some  practical  physical  configura¬ 
tions.  Although  a  relatively  small  amount  of  numerical  data  was  gathered  be¬ 
cause  of  time  limitations,  some  informative  preliminary  results  were  obtained. 

Of  particular  interest  is  a  comparison  that  is  made  with  experimental  data  for 
a  flat  panel  with  free  side  edges,  front  edge  clamped  and  rear  edge  effectively 
pinned  [6j*.  The  experimental  data  give  the  minimum  panel  thickness  required 
to  prevent  flutter  as  a  function  of  Mach  number  (M).  First  mode  flutter  was 
observed  in  all  instances.  For  M  =  1.2  the  present  analysis  predicts  the 
onset  of  first  mode  flutter  in  good  agreement  with  the  experiments.  For 
M  =  1.3  there  is  again  good  agreement  between  theory  and  experiment  on  the 
onset  of  first  mode  flutter.  However,  for  this  latter  Mach  number  the  theory 
shows  that  thicker  panels  will  flutter  in  the  second  mode.  In  other  words, 
the  analysis  indicates  that  first  mode  flutter  is  not  the  critical  condition, 
contrary  to  the  experimental  findings.  These  results  lead  to  the  conclusion 
that  correlation  between  theory  and  experiment  needs  to  be  further  investigated. 


*  Numbers  in  brackets  refer  to  the  bibliography. 
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II.  EQUATIONS  OF  MOTION 


In  this  and  the  following  sections  flutter  equations  of  motion  are 
derived  for  panel  arrays  of  the  type  shown  in  Fig.  1.*  The  panels  are  geo¬ 
metrically  similar  with  length  a  and  width  h  .  The  array  extends  to  an 
arbitrary  number  of  panels  in  the  chordwise  and  spanwise  directions,  and  is 
assumed  to  be  bordered  by  an  inflexible  surface  extending  to  infinity.  The 
upper  surface  (z  >  0)  is  exposed  to  uniform  supersonic  flow  in  the  positive 
x  direction  while  a  uniform,  steady,  pressure  equal  to  the  static  pressure  in 
the  undisturbed  stream  acts  on  the  lower  surface  (acoustic  effects  on  the 
lower  surface  and  membrane  stresses  due  to  static  pressure  difference  are  not 
admitted).  In  the  present  section,  equations  of  motion  are  derived  in  general 
terms.  Specification  of  boundary  conditions  and  exact  formulation  of  the 
aerodynamic  forces  are  presented  in  later  sections. 

Appropriate  equations  of  motion  are  obtained  by  extending  the  derive 
tion  in  [l]  to  include  the  case  of  a  finite  number  of  panels  in  the  spanwise 
direction.  From  small  deflection  plate  theory  the  equation  of  motion  for  the 
array  is 

Djw'(4X  ^  +  2w'(2x''2y'  ^  +  w'^4y  0  +  pshw'*'2t^  +p'(x’,y',t)  =  0  (l)** 


In  Eq.  (l)  w'  is  the  transverse  displacement  in  the  z  direction,  D  the 
flexural  rigidity  of  the  plate,  ps  the  material  density,  h  the  plate 
thickness,  and  p  the  aerodynamic  pressure  excess  on  the  side  z  >  0  . 

Introducing  dimensionless  quantities  x,  y,  w,  p,  and  s  according  to 
x  =  x 1 /a  ;  y  =  y'/b 

w  =  w'/a  i  P  =  P'/pU2  (2) 

s  =  a/b 


Eq.  (l)  in  dimensionless  form  becomes 


*  Figures  and  tables  are  shown  in  Appendix  A. 

**  Differentiation  is  indicated  by  superscripts  in  parentheses  or  brackets. 

This  unconventional  notation  is  adopted  to  help  clarify  subsequent  formu¬ 
lations  . 
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v(4x)  +  2s2v(2x,2y)  +  s^(4^  +  (psha4/D)w(2t ) 

+  (pU2a3/D)p(x,y,t)  =  0  (3) 

Solution  of  (3)  is  obtained  using  the  Ritz-Galerkin  method.  In  the  determina¬ 
tion  of  flutter  boundaries  a  harmonic  solution  is  sought  so  the  deflection  is 
approximated  by 


w  =  elcct  X 
m 


and  the  pressure  is  written 


(4) 


p(x,y,t)  =  elu>tp(x,y-) 


(5) 


where  the  functions  $m(x)  and  ('n(y)  satisfy  the  boundary  conditions  on 
the  panel.*  Following  the  Ritz-Galerkin  method,  (4)  and  (5)  are  substituted 
into  (3).  The  resulting  equation  is  multiplied  through  by  §-(x)tyn(y)  and 
integrated  across  the  length  and  width  of  the  panel  array  to  give, 


z 


^n 


(4x)*_ 


§a  dx 


rN  2  2  f 

<  dy  +  2s2 
Jo  Jo 


$(2x)§_  dx 


*m 


,'.(2y) 


dy 


4  2 
Psha  oq 


^m^m 


P(x,y)yn  dy  dx 


0 


(6) 


Now  define  dimensionless  parameters 


V-  =  Tp  s/p 

(7) 

Z  =  T3E/q(l-v2) 

(8) 

*  The  reason  for  representing  the  spanwise  dependence  of  w  by  a  single 
function  ^n(y)  is  discussed  later. 
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where 


where 


k  = 


coa 

U 


Prior  to  numerical  evaluation  of  (10)  two  subsidiary  steps  are  necessary 
First  the  deflection  functions  $m(x)  and  ^ n (y)  must  be  formulated, 
the  aerodynamic  integral  must  be  put  in  a  form  suitable  for  computation. 
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Then 


III.  THE  DEFLECTION  FUNCTIONS  4m(x)  and  i|[n(y) 


The  deflection  functions  §m(x)  an(l  ^n(y )  must  satisfy  the 
boundary  conditions  on  the  panel  array.  In  the  present  study  two  general 
physical  configurations  are  to  be  analyzed,  each  with  its  own  set  of  boundary 
conditions.  These  configurations  are: 

A.  A  panel  array  as  pictured  in  Fig.  1,  in  which  the  panel  edges  are 

supported  by  elastic  ribs  and  stringers  in  the  chordwise  and 
spanwise  directions, respectively. 

B.  A  panel  array  with  one  panel  in  the  spanwise  direction  (N  =  l), 

with  the  side  edges  free  and  the  spanwise  edges  restrained  by 
stringers  as  in  A. 

The  deflection  functions  are  formulated  for  configuration  A  first. 

Configuration  A  -  Boundary  conditions  for  the  panel  array  supported 
by  elastic  ribs  and  stringers  are  complex  and  the  exact  conditions  are  not 
treated  here.  Instead,  two  simplifying  assumptions  are  introduced  as  follows: 

1.  The  supporting  structure  is  infinitely  rigid  in  bending,  and 

2.  Torsion  is  transmitted  along  the  line  of  a  rib  (stringer)  in 

the  plate  only,  i.e.,  the  rib  (stringer)  offers  resistance 
to  rotation  proportional  to  local  rotation  in  the  plate. 

These  assumptions  retain  most  of  the  important  physical  features  and,  at  the 
same  time,  reduce  the  complexity  of  the  problem.  As  a  consequence  of  the 
assumptions  the  transverse  deflections  at  the  edges  of  the  panels  are  zero 
and  the  bending  moments  are  proportional  to  the  slope  perpendicular  to  the 
edges . 

In  formulating  the  boundary  conditions  along  the  stringers  it  is 
convenient  to  introduce  a  local  coordinate  (Fig.  2). 


9  =  x-(i-l)  ;  jfc-1  £  x  £  & 


and  a  local  chordwise  deflection 

(e) 

5  =  ft  y 

vm  *m 


(id 
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l  =  1/2,3 


L 


(12) 


The  elastic  restraints  against  rotation  at  the  leading  and  trailing 
edge  stringers  are  taken  as  one-half  the  restraint  at  the  intermediate  ones 
to  simplify  the  formulation.  The  boundary  and  compatibility  conditions  follow 
from  the  requirement  of  zero  deflection  and  continuity  of  slope  and  moment  at 
the  panel  edges. 


5m,I<0)  ’ 

Vj(1)  = 0 

I  =  1,2,3  .. 

. .  L  (13a) 

,(«)  - 

I  =  2,3  .... 

L  (13b) 

-  ^x49](°) 

Jo  —  2,3  .... 

L  (13c) 

(29) 

*m,l  <°>  “ 

(9) 

ex*m,l(°) 

(13d) 

(29) 

*m,L  " 

(9),  . 

"  ex^m,l/^ 

(13e) 

where  e  is  a  dimensionless  coefficient  expressing  the  stringer  restraint 
against  rotation.*  Equation  (13a)  requires  zero  deflection  at  the  stringers 
and  (13b)  requires  continuity  of  slope  across  the  stringers.  Equation  (13c) 
is  a  moment  balance  at  the  intermediate  stringers  while  (13d)  and  (I3e)  are 
moment  balances  at  the  bounding  ones. 

Eoundary  conditions  along  the  ribs  are  formulated  in  similar  manner 
by  introducing  local  coordinates  and  spanwise  deflection 


Tj  =  y-(lc-l)  ;  k-1  s  ysk  (14) 

and 


*  In  terms  of  dimensional  quantities,  e„  is  given  by 

X 

T  a 

0  =  _  - 

X  Ddw/dx 

where  Tx  is  the  torque  per  unit  length  in  the  stringer  and  £w/5x  is 
the  slope  perpendicular  to  the  stringer. 
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lc  =  1,2,3  ....  N 


(15) 


The  boundary  and  compatibility  conditions  become 


V<0>  -  V(1>  - 0 


;  k  =  1,2,3 - N  (16a) 


fSL'1)  ■  *$(o> 


j  k  =  2,3  ....  IT  (16b) 


+  «y«SL(1>  -  d)(0)  -  V<£k<0>  -•k-2'3---  "  (16c) 


a’to)  - 


(I6d) 


C>(1)  ’  -  V$<1) 


(16e) 


where  ey  is  a  coefficient  expressing  the  rib  restraint  against  rotation.* 

The  deflection  functions  §m  and  are  taken  to  be  the  natural 

vibration  mode  shapes  of  beams  having  the  boundary  conditions  (13)  and  (16), 
respectively.  The  natural  frequencies  and  mode  shapes  of  a  continuous  beam 
simply  supported  at  equal  intervals  are  given  by  Miles  [ 2 J  .  Miles ’  analysis 
was  extended  to  include  restraint  against  rotation  in  [lj  .  This  development 
in  Ql^  is  repeated  in  Appendix  B  in  slightly  different  form.  It  is  shown  that 
the  deflection  functions  can  be  written  in  the  form 


*  In  terms  of  dimensional  quantities,  € is  given  by 

Cy  =  Tyb/D(3w/Sy) 

where  Ty  is  the  torque  per  unit  length  in  the  rib  and  Bw/dy  is  the 
slope  perpendicular  to  the  rib. 
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(17) 


and 


where 


Cm,Ifm(0) 

+  Dra,  J£’m(1'0 ) 

^n,k 

VnW 

+ 

f  (9)  = 
nr  ' 

sin  Vme  - 

sin 

Vm  sinh  Ym9 

sinh  Ym 

sin  YrJ!  - 

sin 

Yn  sinh  Yn'H 

syn)  - 

sinh  Yn 

(18) 


(19) 


(20) 


Complete  formulations  for  *  r  and  <r  v  are  given  in  Appendix  B. 

Configuration  B  -  Recall  that  this  configuration  consists  of  a 
panel  array  as  shown  in  Fig.  1  having  one  panel  in  the  spanwise  direction 
(N  =  1);  with  the  side  edges  free  and  the  spanwise  edges  restrained  hy 
stringers.  It  is  assumed  that  the  panel  deflects  two -dimens ionally,  i.e., 

Sw/dy  =  0  ,  and  therefore  that  =  1  .  The  chordwise  deflection  function 

§m  is  the  same  as  for  configuration  A. 

With  the  deflection  functions  §m  and  defined,  the  integrals 

of  (9)  can  be  readily  evaluated  in  closed  form.  For  these  evaluations  the 
reader  is  referred  to  Appendix  D. 

One  problem  remains  in  the  development  of  usable  flutter  equations 
from  (10),  namely,  integration  of  the  aerodynamic  pressure  term.  This  in¬ 
tegral  is  evaluated  in  the  following  section. 
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IV,  AERODYNAMIC  INTEGRAL 


It  is  required  to  evaluate  the  integral 


from  (10),  where  p(x,y)  is  the  aerodynamic  pressure  acting  on  the  top  sur¬ 
face  of  the  panel  array.  Since  the  low  supersonic  Mach  number  region  is  of 
particular  interest, the  pressure  is  obtained  from  linearized,  exact,  three- 
dimensional  aerodynamic  theory.  Use  is  made  of  a  result  from  the  analysis 
of  Luke  and  St.  John  [3]  which  states  that  the  perturbation  pressure  on  the 
upper  surface  for  harmonic  motion,  arbitrary  chordwise  deflection  §m(x)  , 
and  sinusoidal  spanwise  deflection  sin  utt0/B  can  be  written  in  the  dimension¬ 
less  form* 

pu  =  i  sin  2Z£  {  4X)  +  jk[(M2-2)/p2l$m  +  j  Pu(S)§m(x-S)  dd  (21) 

'331  J  o  J 


where 

pu(§)  =  e~mi  j-  [(ru2/2)  +  (k2/p4)]jo(ru§)  +  (ru2/2)J2(ru§) 

+  j(2kru/p2)Jx(ru§)  ■  (22) 

K  =  kM/p2  ;  ru2  =  K2  +  (uns/Bp)2  ;  p2  =  M2  -  1 

and  Jn  are  Bessel  functions  of  the  first  kind. 

The  relative  simplicity  of  this  result  comes  from  having  a  deflection 
shape  which  extends  indefinitely  in  the  spanwise  direction.  This  advantage 
can  be  extended  to  panel  arrays  of  finite  spanwise  extent  by  expanding  the 
spanwise  deflection  shapes  in  a  sine  series.  The  total  pressure  can  then 

*  Bb  is  the  dimensional  half  wavelength  of  the  spanwise  deflection  Shape. 
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be  expressed  as  a  superposition  of  terms  similar  to  Pu^m 
let  the  spanwise  deflection  be  represented  by- 


in  (21).  Therefore 


=  Z  Bn,u  sin 
u 


(23) 


The  expansion  (23),  being  an  odd  periodic  function,  gives  flanked  by  ^ 

periodic  reflections.  Therefore  B  is  to  be  chosen  such  that  the  reflections 
are  uncoupled  aerodynamically,  in  order  to  achieve  the  effect  of  an  isolated 
finite  panel  array.*  The  variable  j)  is  the  shifted  coordinate 


0  =  y+X 


where  the  constant  X  depends  on  B.  Criteria  for  the  selection  of  B  and 
X  are  derived  in  Appendix  E. 

From  (21)  and  (23)  the  pressure  on  the  mn'th  deflection  shape  is 


m,n 


^m,n  ]L 
u 


B  P 
n,u5u,m 


(24) 


and  the  general  pressure  tern  becomes 


»<*>*>  ■  z  -  z 

m  m 


qm,n  Z  Bn,upu;m 
u 


(25) 


After  some  manipulation  the  aerodynamic  integral  from  (10)  may  be  written 

J  J  P(x,0)$fi(xHn(0)  d0  dx  =  i  Y_  %,n  |  ^m3 

+  [(m2-2)/P2J  +  Z  Bn,uF(u)Im,m,u]  (26> 


*  For  further  discussion  of  this  point  see  Appendix  E  and  Figs.  3  and  4. 
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where 


R- 

m.m 


$™(x)§-(x)  dx 
*mx  1  m 


(27) 


I_ 

m,m,u 


Pu(§)5m(x-5)^fi(x)  dg  dx 


(28) 


PX+N 

F(u)  =  /  *  sin  HUE.  d0 

4  B 


(29) 


The  integral  R-  is  obtained  in  closed  form  by  mahing  use  of  (17)  and  (19) 
(see  Appendix  vf’.  The  integral  F(u)  can  also  be  readily  obtained  in  closed 
form  hut  it  is  convenient  to  momentarily  postpone  further  discussion  of  this 
integral. 


Two  procedures  can  be  used  to  evaluate  I-  • 

in  y  ul  y  1-1 

1.  When  r  5  and  m  are  small  the  Bessel  functions  in  P  (?)  can 
u  r  t  ^ 

be  approximated  by  a  sum  of  circular  functions  [3 j  such  that 


where 

ar  =  -  (l/q)[(k2/f34)  +  \r2] 

b  =  2k\  /q£2  (31) 

r  r' 

Xr  =  Pu  cos  [(2r-l)/4qjrr 


Then  if  5  (x)  is  expanded  in  a  Fourier  sine  series 
m 
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(32) 


5raW  =  E  \t  sin  (WDx 


j_  can  be  evaluated  in  closed  form.  Details  of  this  evaluation  are  given 

in' Appendix  C.  Roughly  speaking  this  procedure  is  satisfactory  for  m  £  4  , 

L  <:  2  ,  and  I*  5  £  10  . 

2.  When  the  first  procedure  is  not  admissible  it  is  convenient  to 
first  change  the  order  of  integration  to  give 


*m,m,u 


Pu(§) 


§n(x-S)$ni(x)  dX 


(33) 


The  integration  over  x  can  be  completed  in  closed  form.  For  this  purpose 
introduce  another  local  coordinate  similar  to  (ll) 


then 


where 


cp  =  5-U-1)  ;  a-i  £  §  ^  a 


l  n 


I-  =  E  I  P,U-1+?)G-  m  *? 

m,m,u  ^  Jq  n  m,m,X, 


L-i  r  i 

+  E  Pu(A-1+cP)Hm,tM(c?)  d? 

jj=l  Jo 


,*<»>-£  f  * 

A=i  Jep 


-  t  dS 

£=4+1  J  o 


(34) 


(35) 


(36) 
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G-  and  ^  are  integrated  in  closed  form  in  Appendix  C.  The  integra- 

tions^in  (35)  ar e ; evaluat ed  numerically.  The  first  procedure  for  evaluating 
Igj  m  u  is  used  in  a  previously  developed  computer  program  for  calculating 
flutter  boundaries  for  panel  arrays  of  infinite  spanwise  extent.  The  limita¬ 
tions  imposed  by  the  Bessel  function  approximation  (Eq.  (30))  prohibit  its 
use  on  finite  span  arrays.  Consequently,  the  program  developed  here  makes  use 
of  the  second  procedure. 

It  remains  to  evaluate  the  Fourier  coefficients  Bn  u  and  the 
function  F(u)  introduced  in  (26)  and  (29).  These  quantities  depend  on  both 
geometrical  and  aerodynamic  details  and  it  is  thus  appropriate  at  this  point 
to  describe  the  specific  physical  situations  that  are  to  be  analyzed.  There 
are  five  cases  in  all: 

1.  A  panel  array  as  shown  in  Fig.  1  with  a  finite  number  of  panels 
in  both  the  chordwise  and  spanwise  direction.  The  array  is  bordered  by  an 
inflexible  surface  extending  to  infinity  in  all  directions. 

2.  The  same  as  (l)  except  that  the  array  extends  to  infinity  in  the 
spanwise  direction. 

3.  The  same  as  (l)  except  that  the  array  is  flanked  on  the  sides  by 
vertical  walls  (a  wind  tunnel  installation  for  example)  an  arbitrary  distance 
from  the  edge  of  the  array.* 

4.  The  same  as  (l)  except  that  the  array  has  one  panel  in  the  span- 

wise  direction  (N  =  l)  and  the  side  edges  are  free. 

5.  The  same  as  (4)  except  that  the  array  is  flanked  by  vertical 

walls  as  described  in  (3). 

Cases  (l)  -  (3)  conform  to  the  general  geometrical  configuration  A  described 
earlier,  while  (4)  and  (5)  conform  to  configuration  B.  The  effect  of  the 
vertical  walls  in  (3)  and  (5)  is  obtained  by  introducing  image  configurations 
on  each  side  (Fig.  4).  For  detailed  evaluation  of  Bn^u  and  F(u)  the 
reader  is  referred  to  Appendix  E. 


*  The  analysis  is  completed  for  this  option  but  the  option  was  not  programmed 
because  of  time  limitations. 
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V.  FLUTTER  EQUATIONS  AND  METHOD  OF  SOLUTION 


Flutter  equations  are  obtained  from  (10)  and  (26) 

Z!  <  ! 7  [YmJm,ms  +  2s2KE,mT  + 

m  L^4 


+  21  Bn^uF(u)Im^m^u  f  ~  ^ 


In  (37)  the  stiffness  terms  are  multiplied  by  (l+jg)  to  account  for  struc¬ 
tural  damping. 

It  is  convenient  to  write  (37)  in  matrix  form.  Therefore >  let 


Vm^m^m  +  ^s^^m(,m^/^  +  s 


Cm,m|  =  |  Rm,m  +  jJm,mk(M2-2)/p2  +  Am,m  +  j^,mj 
'  =  Real  <  Bn^uF(u)lm^m^u/s^ 

=  Imag  Bn^uF(u)ljn;injU/S ^ 


Then  (37)  can  be  written 


where 


[‘Imj  1 


(l+jg)  E^m  -  p,k2J; 


is  the  flutter  vector. 


i,m  +  |  Cm, ml 


Im  p  ~  0 
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Equation  (39)  represents  a  set  of  simultaneous,  homogeneous,  algebraic 
equations.  For  a  nontrivial  solution  it  is  necessary  to  have 


Det  jfj  %,m  “  +  |  =  0  (40) 


Equation  (40)  is  a  concise,  mathematical  statement  of  the  flutter  problem. 

The  computer  program  must  first  calculate  the  elements  of  the  flutter  matrix 
in  terms  of  selected  ’'free"  parameters,  and  then  find  what  combinations  of 
free  parameters  satisfy  (40).  Setting  up  the  matrix  is  straightforward  and 
presents  no  difficulty.  The  process  of  satisfying  (40),  however,  requires 
some  discussion. 

On  the  surface  it  appears  that  the  free  parameters  might  be  any  of 
[i  >  Z?  k  and  g  •  The  reduced  frequency,  k  ,  however,  is  inadmissible  since 
it  is  contained  within  Cfn,m  in  complicated  transcendental  form.  Further,  it 
is  desirable  to  make  g  an  input  parameter,  i.e.,  to  be  able  to  perform 
flutter  calculations  for  specific  values  of  g  .  This  leaves  only  the  two 
real  parameters  p,  and  Z  ,  and  since  the  flutter  determinant  of  (40)  is 
complex  it  is  necessary  that  both  p,  and  Z  be  retained  as  free  parameters. 
The  flutter  problem  is  therefore  reduced  to  one  of  finding  combinations  of 
real  p,  and  Z  which  satisfy  (40),  for  given  input  values  of  g  and  k  . 

There  is  no  convenient  closed-form  type  of  solution  to  the  posed 
problem.  A  conventional  approach  is  to  insert  specific  real  values  for,  say 
[i  ,  and  solve  for  the  Z's  which  satisfy  (40).  The  Z's  so  obtained  are 
complex  in  general  and  it  is  necessary  to  interpolate  on  p,  to  find  where 
the  imaginary  part  of  Z  vanishes.  The  process  is  complicated  by  the  fact 
that  for  each  p,  there  are  m  values  of  Z  which  satisfy  (40).  This  compli¬ 
cation  ultimately  leads  to  the  necessity  to  make  successive  computer  runs  to 
obtain  flutter  points  for  one  set  of  g,k  values. 

In  the  present  analysis  a  solution  technique  is  used  which  reduces 
computation  time  by  obtaining  flutter  points  in  a  single  computer  run.  The 
method  proceeds  as  follows.  For  a  given  value  of  p,  the  characteristic 
equation  of  the  determinant  is  obtained.  This  equation  is  a  polynomial  in  Z 
with  complex  coefficients.  Assuming  Z  to  be  real,  the  real  and  imaginary 
parts  are  separated  into  two  polynomials 


a^Zm  + 


a^Z  +  a0  =  0 


(41) 


+ . bxz  +  b0  =  0 

Any  value  of  Z  (with  its  associated  p,  )  which  satisfies  both  of  Eqs.  (41) 
constitutes  a  flutter  point,  A  modified  flutter  criterion  is  therefore  ob¬ 
tained,  namely  that  Eqs.  (41)  have  at  least  one  common  root,  Now  (41)  will 
have  a  common  root  if  and  only  if  the  Sylvester  determinant  formed  from  the  aj_ 
and  bj_  is  zero  [lj .  This  Sylvester  determinant  is  in  essence  a  single-valued 
function  of  [x  in  contrast  with  the  multi-valued  function  that  arises  in  the 
conventional  approach  to  the  solution.  In  the  present  computer  program  the 
Sylvester  determinant  is  calculated  for  a  specified  array  of  [x  values,  and 
an  interpolation  procedure  is  then  used  to  find  the  exact  values  of  [x  for 
which  it  vanishes.  These  jjl  ,  with  associated  Z  values,  constitute  flutter 
points.  The  complete  process  is  accomplished  in  a  single  computer  run. 

A  computer  program  has  been  developed  which  uses  the  solution  tech¬ 
nique  just  described.  Numerical  considerations  dictate  the  size  of  flutter 
system  that  can  be  handled  by  the  program,  i.e.,  the  number  of  chordwise  modes, 
m  ,  that  can  be  used.  The  present  program  operates  satisfactorily  in  general 
for  m  up  to  six.  Specially  tailored  routines  are  used  for  scaling,  determi¬ 
nant  evaluation  and  characteristic  equation  expansion.  It  is  anticipated  that 
further  refinements  in  these  areas,  accomplished  through  selective  application 
of  known  computational  techniques,  would  extend  the  range  of  general  applica¬ 
tion  to  a  greater  number  of  chordwise  modes. 
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VI.  NUMERICAL  RESUITS 


The  computer  program  obtains  flutter  points  in  terns  of  the  two 
general  parameters 


and 


l/u,  -  -  air-panel  mass  ratio  , 

hps 


jl/3  _  h  E 

a  _q(l-u2) j 


,1/3 


=  stiffness -dynamic  pressure  ratio 


Flutter  boundaries  can  be  constructed  in  the  l/p,  -  Z  '  plane  by  making 
successive  calculations  for  judiciously  selected  values  of  reduced  flutter 
frequency  k  •  When  only  the  critical  boundary  is  sought  it  is  still  important 
that  a  range  of  k  be  examined  which  is  extensive  enough  to  locate  all  the 
flutter  boundaries  in  the  l/p,  -  Z1'3  plane.  If  this  is  not  done  the  bound¬ 
aries  identified  as  critical  may  be  open  to  question. 


Most  of  the  effort  during  this  contract  period  went  into  development 
of  the  computer  program.  Consequently,  a  limited  quantity  of  numerical  data 
was  obtained.  Critical  flutter  boundaries  were  computed  for  panel  arrays  with 
one  chordwise  bay  (Figs,  6  and  7).  Calculations  were  also  carried  out  for 
panel  arrays  with  three  chordwise  bays.  These  calculations  were  exploratory 
ones  designed  to  test  the  adequacy  of  the  computer  program.  Finally,  calcula¬ 
tions  were  made  for  panels  with  free  side  edges,  some  in  the  presence  of  out¬ 
board  vertical  walls,  for  comparison  with  experimental  data. 


Flutter  Boundaries  for  an  Array  of  Panels 


Critical  flutter  boundaries  obtained  in  a  four  mode  analysis  for  a 
panel  array  consisting  of  one  chordwise  bay  and  an  infinite  number  of  spanwise 
bays  are  shown  in  Figs.  6  and  7.  Three  edge  conditions  are  represented  at 
each  of  three  Mach  numbers,  for  an  aspect  ratio  of  4  and  structural  damping 
coefficient  of  0.01.  The  edge  condition  c  =  0  is  the  pinned  condition  while 
e  =  1,000  closely  approximates  the  clamped  case.* 


*  The  pinned  edge  results  were  obtained  earlier  under  this  contract  and  are 
included  in  [jbJ.  They  are  repeated  here  for  comparison  purposes. 
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For  M  =  1.25  and  2  the  flutter  vectors  associated  with  the  boundaries 
are  dominated  by  the  first  vacuum  mode.  Coupling  with  the  other  modes  is  pres¬ 
ent  in  both  cases  but  is  more  pronounced  for  M  =  2  .  For  M  =  /*2~  (Fig.  7) 
the  boundaries  for  e  *  10  and  1,000  are  also  first  mode  dominant.  For  €  =  0 
the  critical  boundary  is  composed  of  three  distinct  regions.  The  portion  to 
the  right  is  third  mode  dominant,  the  short  loop  is  fourth  mode  dominant,  and 
the  remainder  is  first  mode  dominant. 

To  interpret  the  results  in  Figs.  6  and  7  it  is  best  to  view  values 
on  the  abscissa  ( l/p, )  as  denoting  specific  combinations  of  panel  material, 
panel  thickness,  and  altitude.  Then  the  ordinate  values  (Z1/3)  of  points  on 
the  boundaries  are  inversely  proportional  to  the  minimum  free  stream  dynamic 
pressure  that  will  cause  flutter.  Thus  larger  Z1/3  values  correspond  to  less 
stable  conditions.  It  is  seen  in  Figs.  6  and  7,  that  in  all  cases,  decreasing 
c  has  a  destabilizing  effect,  i.e.,  flutter  occurs  for  larger  values  of  zV  . 
It  is  also  apparent  that  the  effect  of  varying  g  diminishes  with  increasing 
Mach  number.  Finally,  comparison  of  the  figures  shows  that  increasing  the 
Mach  number  has  a  stabilizing  effect. 

Six  mode  analyses  were  conducted  for  pinned  edge  panel  arrays  with 
three  chordwise  bays  and  an  infinite  number  of  spanwise  bays.  These  calcula¬ 
tions  were  intended  primarily  to  assess  the  adequacy  of  the  computer  program. 
They  were  not  extensive  enough  to  permit  the  identification  of  critical 
boundaries  and  therefore  no  detailed  results  are  given  here. 


Flutter  Boundaries  for  a  Single  Finite  Panel 

Some  calculations  were  carried  out  for  comparison  with  the  experi¬ 
mental  data  in  [e]  and  the  analytical  results  in  [7].  The  experimental  data 
are  for  a  panel  with  free  side  edges,  front  edge  clamped  and  rear  edge  effec¬ 
tively  pinned.  There  is  a  gap  between  the  side  edges  and  the  wind  tunnel  wall 
equal  to  6.2  per  cent  of  the  panel  span.  Aspect  ratios  vary  between  0.95  and 
1.06.  The  analysis  £7]  treats  the  case  of  a  panel  with  side  edges  free  and 
does  not  include  the  effect  of  the  wind  tunnel  wall.  In  both  [6]  and  [V], 
results  are  presented  in  the  form  of  thickness-to-pre vent -flutter  as  a  function 
of  Mach  number. 

Two  types  of  calculations  were  made  here,  namely,  for 

(a)  A  panel  with  side  edges  free  surrounded  by  a  quiescent 
surface  as  in  QJ ,  and 

(b)  A  panel  with  side  edges  free,  flanked  by  vertical  sidewalls 
to  simulate  a  wind  tunnel  installation. 
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Four  mode  analyses  were  conducted  using  the  first  four  natural  vibration  modes 
of  a  two-dimensional  beam.  Twenty  terms  are  carried  in  the  Fourier  expansion 
of  the  spanwise  deflection  shape.  (Typical  twenty  term  representations  are 
shown  in  Fig.  8).  All  calculations  are  for  an  aspect  ratio  of  unity  and  a 
structural  damping  coefficient  of  0.01  (approximately  the  mean  experimental 
values).  The  experimental  edge  conditions  are  not  duplicated  exactly  since 
the  present  analysis  assumes  that  edge  conditions  at  the  front  and  rear  of 
the  panel  are  the  same.  Instead,  results  are  obtained  for  pinned  (e  =  0)  and 
clamped  (e  =  1,000)  conditions  to  bracket  the  experimental  conditions. 

The  calculations  were  carried  out  first  for  configuration  (a)  at 
M  =  1.2  .  The  following  critical  thickness  ratios  (Tcrj_t)  an(^  frequencies 
were  obtained. 

Dominant 


e 

Tcrit 

k* 

Mode 

0 

0.00465 

0.134 

1 

1,000 

0.00424 

0.275 

1 

These  points  agree 

remarkably  well  with  the 

experimental  results  for  M  =  1.2  , 

namely, 


Tcrit  "**  0*00438 
k  =  0.187 


For  this  same  case  the  analysis  b]  gives  Tcrj_t  ”  at  a  frequency 

"slightly  below  the  second  natural  frequency. "  It  is  seen  that  the  present 
results  agree  better  with  the  experimental  data  than  with  £7],  even  though 
the  tunnel  walls  are  not  simulated. 

Calculations  were  next  carried  out  for  configuration  (a)  at  M  =  1.3  . 
The  following  results  were  obtained. 


*  Exact  values  of  k  are  not  determined  by  calculation  since  Tcrj_^  is 

obtained  by  finding  where  the  critical  boundary  intersects  the  line  in 
the  i/jjl  -  z1'3  plane  denoting  experimental  conditions.  The  values 
given  for  k  come  from  linear  interpolation  between  calculated  points 
on  each  side  of  the  intersection. 
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Dominant 

Tcrit  k  Mode 


0.00564  0.61  2 

0.00426  0.73  2 


For  this  case  the  experimental  results  are 

Tcrit  =  °*00335 
k  =  0.165 

and  the  analysis  [7  3  gives  Tcrit  =  0.0038  and  frequency  "slightly  below  the 
third  natural  frequency. " 

The  failure  to  agree  with  the  experimental  data  for  this  case 
prompted  a  closer  look  at  the  numerical  data.  It  was  noted  that  subcritical 
boundaries  associated  with  the  first  chordwise  mode  are  present  in  the  region 
of  general  interest.  It  was  not  possible  to  run  additional  calculations  to 
determine  exactly  where  these  boundaries  intersect  the  experimental  condition 
line.  However,  reasonable  extrapolation  yields  the  following  thickness  ratios 


€ 

0 

1,000 


€  T 

0  0.0036 

1,000  0.0032 


The  corresponding  reduced  flutter  frequencies  cannot  be  obtained  by  extrapola¬ 
tion  but  it  can  be  seen  that  they  are  less  than  0.25.  These  values  of  t  and 
k  again  compare  well  with  the  experimental  data. 

Calculations  were  next  carried  out  for  configuration  (b)  to  assess 
the  effect  of  the  wind  tunnel  walls.  The  following  case  was  investigated. 


M  =  1.3 


aspect  ratio  =,1 

g  =  0.01 
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gap  to  vail  =  5  per  cent  of  span 
edge  conditions  =  pinned;  clamped 

The  5  per  cent  gap  is  smaller  than  the  actual  gap  in  the  experimental  setup 
(6.2  per  cent).  It  was  so  taken  to  allow  for  the  boundary  layer  displacement 
thickness  effect.  The  folloving  results  were  obtained 

Dominant 

G  Tcrit  k  Mode 


0  0.0060  0.65  2 

1,000  0.0046  0.80  2 

The  presence  of  the  vails  is  seen  to  increase  Tcrit  by  about  10  per  cent, 
i.e.,  the  vails  have  a  slight  destabilizing  effect.  Examination  of  the 
numerical  data  again  revealed  the  presence  of  subcritical  boundaries  associ¬ 
ated  with  the  first  mode,  but  extrapolation  of  the  boundaries  vas  not  attempted. 

The  results  that  have  just  been  described  are  summarized  belov  for 

clarity. 


Source  of  Mach 

Edge* 

Thickness 

Dominant 

Gap  to* 

Critical 

Data 

No. 

Cond. 

Ratio 

Freq.. 

Mode 

Wall 

Cond. 

Present 

1.2 

P 

0.00465 

0.134 

1 

Yes 

ii 

1.2 

C 

0.00424 

0.275 

1 

- 

t> 

[6] 

1.2 

P-C 

0.00438 

0.187 

1 

6.256 

n 

W 

1.2 

P-C 

0.00530 

- 

2 

- 

n 

Present 

1.3 

p 

0.00564 

0.61 

2 

- 

n 

it 

1.3 

c 

0.00426 

0.73 

2 

- 

n 

[6 ] 

1.3 

P-C 

0.00^35 

0.165 

1 

6.256 

ti 

[7] 

1.3 

P-C 

0.0038 

- 

3 

- 

it 

Present 

1.3 

P 

0.0036 

k<0 . 25 

1 

- 

No 

n 

1.3 

c 

0.0032 

k<0 . 25 

1 

- 

i> 

n 

1.3 

p 

0.0060 

0.65 

2 

556 

Yes 

u 

1.3 

c 

0.0046 

0.80 

2 

55 6 

11 

*  Under 

edge  conditions, 

P  indicates  pinned, 

C  clamped,  and 

P-C  the 

experimental  conditions  of  front  clamped  and  rear  pinned.  The  gap  to  wall 
is  given  in  per  cent  of  panel  span. 
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The  above  results  must  be  considered  as  preliminary,  subject  to  a 
more  definitive  numerical  study.  Calculations  should  be  made  for  additional 
reduced  frequencies  to  insure  that  the  points  identified  as  critical  are  truly 
critical.  Experience  has  shown  that  it  is  sometimes  easy  to  miss  a  complete 
boundary,  which  may  be  the  critical  one,  even  though  a  well  chosen  range  of 
frequencies  is  used.*  In  this  light  the  case  M  =  1.2  should  be  thoroughly 
examined  for  the  presence  of  a  critical  second  mode  boundary.  Additional 
calculations  should  also  be  made  to  examine  the  convergence  of  the  Galerkin 
procedure  and  of  the  Fourier  expansion  of  the  spanwise  deflection  shape. 


Estimates  of  Computing  Time 


During  the  course  of  the  calculations,  estimates  were  made  of  com¬ 
puting  time.  To  present  these  estimates  in  a  meaningful  way,  it  is  necessary 
first  to  describe  the  calculation  procedure.  A  set  of  input  parameters  is 
selected  first.  Included  are 

1.  Geometry  Code** 

2 .  Aspect  Ratio 

3 .  Mach  number 

4.  Edge  conditions  (c) 

5 .  Number  of  chordwise  modes 

6 .  Number  of  terms  in  expansion  of  spanwise  deflection  shape 

7 .  Damping  coefficient,  g 

8 .  Frequency,  k 

A  run  on  the  computer  consists  in  loading  the  data  corresponding  to  inputs  1 
through  7  and  then  cycling  through  a  sequence  of  frequency  values.  The  total 


*  This  difficulty  of  choosing  reduced  frequencies  is  not  unique  to  the  pres¬ 
ent  analysis;  it  is  common  to  all  analyses  that  make  use  of  exact,  three- 
dimensional  aerodynamic  theory. 

**  Geometry  code  identifies  the  situation  to  be  analyzed;  for  instance,  an 
isolated  panel  with  free  side  edges. 
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time  for  a  run  is  comprised  of  two  separable  parts.  There  is  a  setup  time 
during  which  matrices  that  do  not  depend  on  frequency  are  calculated  and 
stored,  and  there  is  the  time  required  to  process  each  frequency. 

The  number  of  flutter  points  that  are  obtained  varies  from  frequency 
to  frequency  and  the  running  time  varies  accordingly.  For  the  four  mode  cal¬ 
culations  that  were  made  for  the  single  panel  with  free  side  edges  the  follow¬ 
ing  averages  were  obtained. 

Set-up  time  -  0.5  min. 

Time  per  frequency  -  0.4  min. 

Flutter  points  per  frequency  -  3 

A  few  runs  made  using  six  chordwise  modes  indicate  that  setup  time  approxi¬ 
mately  doubles  and  running  time  for  each  frequency  increases  by  a  factor  of 
about  2- l/2. 

For  the  runs  using  six  chordwise  modes  for  a  panel  array  with  three 
chordwise  bays,  the  setup  time  was  about  4  min.  and  the  running  time  for  each 
frequency  about  1  to  l-l/4  min.  The  average  number  of  flutter  points  per 
frequency  was  between  4  and  5. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  computer  program  has  been  developed  which  has  three  significant 

features : 

1  It  uses  a  solution  technique  that  obtains  flutter  points  in  a 
single  computer  run. 

2  It  can  handle  up  to  six-by-six  flutter  determinants. 

3  It  can  be  used  to  investigate  a  variety  of  physical  situations, 
including  that  of  a  single  panel  installed  in  a  wind  tunnel. 

Preliminary  numerical  results  suggest  the  possibility  of  obtaining  meaningful 
comparison  with  experiment. 

The  computer  program  should  be  exploited  by  using  it  in  a  definitive 
parametric  investigation  of  the  flutter  of  flat  panel  arrays. 

An  attempt  should  be  made  to  extend  the  capability  of  the  computer 
program  to  ten-by-ten  systems,  by  refining  the  implicit  numerical  operations. 
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APPENDIX  A 

FIGURES  1  THROUGH  8  AND  TABLE  I 
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TABLE  A-I 


COMPARISON  BETWEEN  PARTIALLY  CLAMPED  FREQUENCIES 
AND  CLAMPED  FREQUENCIES  (L  =  l) 


Mode  e 


1  0 

2  0 

3  0 

4  0 


Partially  Clamped  Frequency  (e) 
Clamped  Frequency  (e  =  «) 

0.441 

0.640 

0.735 

0.790 


1 

2 

3 

4 


10 

10 

10 

10 


0.  772 
0.810 
0.838 
0.859 


1 

2 

3 

4 


100 

100 

100 

100 


0.963 

0.964 

0.965 

0.966 


1 

2 

3 

4 


1,000 

0.996 

1,000 

0.996 

1,000 

0.996 

1,000 

0.996 
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- X 

^  e 


l  =  1  l  =  2 


£ 


(a)  Chordwise  Direction 


■**  y 


■i  J- 


k  =  1 


k  =  2 


(b)  Spanwice  Direction 


Fig.  2  -  Local  Coordinates 


N 


(a)  Actual  Deflection  Shape 


(b)  Analytical  Deflection  Shape 


(c)  Plan  View  Showing  Maximum  Spanwise  Excursion  of 
Forward  Facing  Mach  Cone  (Shaded  Area) 


Fig.  3  -  Geometrical  Aspects  of  Spanwise  Deflection  Expansion 
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Planes  of  Flow  Symmetry  Which  Simulate  Vertical  Side  Wall 


Image  Panel  Configuration 


(rl)  SOD 
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Fig.  5  -  Natural  Frequencies  for  a  Uniform  Beam  on  Many  Support 
with  Elastic  Restraint  Against  Rotation 


-Half  Wavelength  Shown) 


APPENDIX  B 

EQUATIONS  OF  MOTION  OF  A  UNIFORM  BEAM  ON  MANY  SUPPORTS  V7ITH 
ELASTIC  RESTRAINT  AGAINST  ROTATION 
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In  this  Appendix,  Miles'  analysis  of  a  continuous  beam  on  equally 
spaced  supports  ,  is  extended  to  include  the  effects  of  elastic  restraints 
against  rotation  at  the  supports.*  The  equation  of  motion  of  a  uniform  beam 
in  harmonic  motion  is 


$(4x)  =  x4* 


(B-l) 


In  the  formulation  of  the  boundary  conditions  for  a  beam  supported  at  x  =  0, 
a,  2a,  ....  La,  it  is  convenient  to  introduce  a  dimensionless  local  co- 
ordinate 


e  =  -  -  (jj-l)  (jj-l)a  £  X  s  ia  (B-2) 

a 


so  that 


5(0)  =$_(0 )  £  =  1,2,3  ....  L 

X/ 


(B-3) 


where  y4  =  X4a4  • 


Taking  the  elastic  restraints  against  rotation  at  the  intermediate  supported 
edges  twice  as  large  as  the  elastic  restraint  at  the  end  points,  the  boundary 
conditions  are 


5j(°)  =  5.(1)  =  0 

£  =  1,2,3  ....  L 

(B-4) 

46>(d  -  4e)(o) 

jfc-i  & 

I  =  2,3,4  ....  L 

(B-5) 

,<“>(!)  *  •.?><!>  -  *<S9><0>  -  elf  >< 0) 
a-i  A-i  x  i 

J  =  2,3,4  ....  L 

(B-6) 

•<*>«»  -  elf  >(0)  ;  42S>(1)  -  • 

-  c*£9)(l) 

(B-7 ) 

The  development  given  here  is  nearly  identical 

with  that  given  in  Appendix 

B  of  [l]  ,  the  principal  difference  being  in  the  general  form  that  is 
selected  for  the  deflection  shapes. 
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where  e  is  a  dimensionless  coefficient  expressing  the  stiffness  of  restraint 
of  the  edge  supports.* 

A  general  solution  of  (B-3)  satisfying  (B-4)  can  he  written  as 

§-(9)  =  C_[sin  Y0  -  (sin  y  sinh  Y0)/sinh  yJ 
&  l 1 

+  D-jsinY(l-Q)  -  ^sin  Y  sinh  Y(l-9  )j/sinh  yJ  (B-8) 


Substitution  of  (B-8)  into  (B-5),  (B-6),  and  (B-7)  gives 


Cj  -  pDjr  =  pC J_1  -  Dj_i  Z  -  2,3,4 - L 

(B-9 ) 

qCj-+  Bj  =  Cj_x  +  qD^  I  ~  2,3,4 -  L 

(B-10) 

qCi  +  0-^=0 

(B-ll) 

CL  +  =  0 

(B-12) 

where 

p  =  jjsinh  y  -  sin  yj  |sinh  y  cos  y  -  sin  y  cosh  yj 

(B-13) 

g  =  Fay  sin  y  sinh  y  -  e(sinh  y  cos  y  -  sin  y  cosh  y) 

]-x 

jj:(sinh  y  -  sin  y) 

]  (B-14) 

It  is  convenient  to  note  for  later  use,  that  the  conditions  q. 

=  tl  give  modes 

that  are  identical  with  those  of  a  single  span  of  length  a  ,  namely 


*  See  the  report  body  for  details  of  e  • 
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$_(0)  =  q^  jsin  yQ  -  (sin  y  sinh  yQ)/sinh  yj 

-  q^  jsin  y(l-0)  -  ^sin  y  sinh  y(l-0)^/sinh  yj  (B-15) 

It  can  be  readily  verified  that  (B-15)  satisfies  (B-9)  through  (B-12). 

Continuing  with  the  general  development,  Eq.  (B-9)  and  (B-10)  can 
be  operated  upon  to  give 


Cjj(l-Pq)  =  2pD^  -  Dj.-^l+pq) 

(B-16) 

D£+l  "  2  cos  |i  Dj  +  Di_1  =  0 

(B-17) 

where 

cos  p.  =  p+<3- 

(B-18) 

pq+i 

The  general  solution  of  (B-17)  is 

D-  =  A  cos  J[i  +  B  sin  J\l  (B-19) 

and  from  (B-16)  there  follows 

C-  =  — —  <2p  A  cos  &\L  +  B  sin  Jp, 

&  1-pq  |  L  J 

-  (l+pq)[A  cos  (2-1  )n  +  B  sin  (I-l)n]j  (B-20) 

Substituting  (B-19)  and  (B-20)  back  into  (B-ll)  and  (B-12)  gives,  after  some 
algebraic  manipulation, 

A(cos  p,-q)  +  B  sin  p,  =  0  (B-2l) 
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(B-22) 


B  sin  L\i 


g+pg)(i-(i2)l,=  o 

1  +  P<1  J 


For  a  nontrivial  solution  it  is  necessary  to  have  sin  Ip.  =  0  which  gives 


A  =  sin  p  (b-23) 

B  a  -  cos  p 


Any  value  of  p  which  satisfies  the  condition  sin  Lp  =  0  can  be  used.  Later 
considerations  of  the  frequency  spectra  show  that  a  convenient  form  is 

1%  ■  <4  +  T1) 

VB  -  Y<%>  (B-25> 


where  m  »  1,2,3  _ ,  and  ym  is  the  solution  of  (B-18)  corresponding  to 

p  .  Now  p  satisfies  sin  Lp^.  =  0  so  from  (B-8),  (B-16),  (B-17),  and 
nn  .  in  LU 

(B-23)  the  modes  are 


§ 


sin  jJPju  -  qm  sin  (jl-l)Pm 


<Jm  "  cos  h 


m 


sin  (j^-l)pm  -  sin 

<3m  "  cos  ^m 


f  (1-9) 
mv  ' 


where 

fffl(0)  =  sin  Ym0  -  (sin  Ym  sinh  ym9 )/sinh  Ym 


(B-26 ) 


(B-27) 


"  4(Ym) 


(B-28) 
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If  is  a  multiple  of  tt  the  terms  sin  fpm  and  sin  (jj-l)^  in 

(B-26)  are  zero.  Equation  (B-18)  has  two  possible  solutions  for  this  case, 
namely 


q.  =  cos  \i  =  -1 

;  p  =  arbitrary 

p  =  COS  [L  -  tl 

;  <q  =  arbitrary 

For  the  first  solution,  q  =  tl,  (B-26)  is  indeterminate.  However,  it  was 
pointed  out  earlier  that  for  q  =  tl  the  modes  are  identical  with  those  of  a 
single  span  of  length  a  .  It  follows  from  (B-18),  (B-24),  and  (B-9)  through 

(B-12)  that  the  cases  m  =  1,  2L+1,  4L+1  -  ,  correspond  to  q  =  -1  so  that 

from  (B-15) 

V,|<0)  *  +  fm<l-S>]  (B-a9) 


and  if  m  =  L+l,  3L+1,  5L+1  .  ,  q  =  1  so  that 


$  -(e)  =  f  (e)  -  f  (i-9) 

m,  V  nr  m 


For  the  second  solution  of  (B-18),  i.e.,  p  =  tq  ,  (B-26)  is  degenerate  and 
this  solution  is  discarded. 

The  free  vibration  frequencies  follow  from  the  Ym  which  in  turn 
are  determined  by  (B-18)  and  (B-24).  The  permissible  values  of  y  and  p  are 
illustrated  in  Fig.  5  for  e  =  0  and  10  .  The  smallest  admissible  value  of  y 
is  obtained  for  cos  p  =  -1  ,  corresponding  to  the  smallest  y  obtained  from 

(B-18)  and  (B-24)  with  m  =  1  .  Setting  m  =  2,3  .  L  yields  successively 

increasing  values  of  y  ,  with  Yl  <  1-5056  tt  (the  lowest  frequency  of  a  clamped - 
clamped  beam).  Setting  m  =  L+l  gives  y  =  1.5056  tt  but  this  corresponds  to 
the  solution  p  =  1  of  (B-18)  which  is  degenerate.  This  solution  is  there¬ 
fore  discarded  as  trivial.  The  value  m  =  L+l  gives  a  second  value  y  > 

1.5056  tt  which  is  a  permissible  value,  corresponding  to  the  solution  q  =  1 
of  (B-18).  Successively  increasing  values  of  y  are  obtained  by  setting 

m  =  L+2,  L+3,  .  2L  .  Setting  m  =  2L+1  yields  a  trivial  solution 

Y  =  2.4998  tt  corresponding  to  p  =  -1  and  a  permissible  solution  y  >  2.4998  tt. 
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Additional  frequencies  are  obtained  as  indicated,  with  the  trivial  or  de¬ 
generate  solutions  being  successively  the  ascending  frequencies  of  a  clamped- 
clamped  beam. 


Figure  5  shows  that  the  natural  frequencies  for  the  beam  are  contained 
within  discrete  frequency  bands.  These  bands  are  independent  of  the  number  of 
bays,  L  ,  and  the  width  of  the  bands  decreases  with  increasing  e  ,  i-e.,  as 
the  clamped  condition  is  approached.  For  a  beam  of  L  bays  there  are  L 
discrete  natural  frequencies  in  each  frequency  band.  Thus,  for  an  infinitely 
extending  beam  the  natural  frequencies  form  a  continuous  spectrum  throughout 
each  band. 


Table  I  (Appendix  A)  gives  some  insight  into  the  physical  significance 
of  e  .  In  the  table  the  first  four  frequencies  for  a  one  bay  beam  (L  =  l) 
with  e  -  0,  10,  100,  and  1,000  are  compared  respectively  to  the  first  four 
frequencies  for  a  clamped -clamped  beam  (e  =  »)  .  It  is  seen  that,  frequency- 
wise,  c  =  10  lies  approximately  midway  between  the  pinned  case  (e  =  0)  and 
the  clamped  case,  whereas  e  =  1,000  very  closely  approximates  the  clamped 
condition. 


To  summarize,  the  chordwise  deflection  functions  used  in  the  flutter 
analysis  are  as  follows:* 


C  -f*  (0 ) 


+  D  -f  (1-6) 


(B-31) 


where 

cm,z =  sin  K  -  %  sin 

=  sin  -  %  sin 

(B-32) 

fje)  =  sin  Yme  -  (sin  ym  sinh  Ym0  )/sinh  ym 

(B-33) 

q .  =  q(ym) 

m  m 

(B-34) 

/ 


-4 


in  (B-26)  is  discarded. 


*  The  multiplying  factor 


1 


with  the  special  cases 


=  DmJ=  (-1) 


SL- 1 


(B-35) 


for  m  =  1,  2L+1,  4L+1, 


and 


Cm,£ 


1 


Dm,je 


-1 


(B-36) 


for  m  =  Ih-1,  3L+1,  5L+1, 


In  all  cases  frequencies  are  obtained  as  described  earlier,  taking  e  =  ex  . 

The  spanwise  deflection  function  (y)  is  not  formulated  here  since 
it  can  be  deduced  from  the  formulations  for  §  (x)  by  straightforward  manipu¬ 
lation  of  the  notation,  taking  e  =  €y  • 
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APPENDIX  C 

EVALUATION  OF  THE  AERODYNAMIC  INTEGRAL  Ig 
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From  Eqs.  (22)  and  (28)  of  the  report 


I_* 

m,m,u 


/  /  Pu^)5m(x_5)$m(x)  d5  dx 


where 

Pu(s)  =  e'JKM§  |-[(ru2/2)  +  (k2/p4)]j0(ru§)  +  (ru2/2)J2(ru§) 

+  j(2kru/p2)J1(ru?) j 

K  =  kM/p2  }  ru2  =  K2  +  (uns/Bp)2  p2  =  M2-l  j2  =  -1 
Two  procedures  for  evaluating  I-  ,,  are  outlined  in  the  report.  The  pur- 

m,  m  y  u. 

pose  of  this  Appendix  is  to  present  these  procedures  in  detail. 


Procedure  1 


This  procedure,  the  more  approximate  of  the  two,  makes  use  of  two 
basic  approximations.  First,  the  Bessel  functions  in  Pu(§)  are  approximated 
by  sums  of  circular  functions  [3j  • 


J0(ru5)  =  -  Z  cos  (Xr§)  +  2J4<i(ru5) 

<1  r=l 


Jl(ru5)  =  i  Z 

*1  r=l 


2r-l 

cos  |  =i-=-  |n 
4q 


sin  (Xr§) 


(C-l) 


["  J4q-l(ru?)  +  J4q+1  <r„!>] 
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where 


x_  =  r„  cos 


(¥) 


rr 


(C-2) 


The  term(s)  appearing  outside  of  the  summation  sign  give  an  estimate  of  the 
error  in  the  representation  and  consequently  are  omitted  in  the  analysis.  In 
the  numerical  work,  ^(Fu? )  replaced  by  means  of  the  identity 

2J  (r e)  ,  , 

J2(ru§)  = . -u—  -  J0(ru5) 


The  second  basic  approximation  introduced  in  this  procedure  is  the 
expansion  of  the  chordwise  deflection  function,  ,  in  a  truncated  sine  series 


=  T  \t sin 


t-rrx 

L 


(C-3) 


Substitution  of  these  approximations  gives 


I-  =  L2  X  2  A  +A-tZ+7 
m,m,u  ~  ^  m,t  m,t  t,t,u 

"G  O 


(C-4) 


where 


-  =  i  E,  6,  t  +  -  <-F.  e.  *?  -  (-1) 

;,t,u  2  t  t,t  L  l  t  t,t 


*  Ht 


4 


(C-5) 


6,  t  =  0  if  t  ^  t 

(0-6) 

=  1  if  t  =  t 
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G  =  0  if  (t+t )L  is  even 
t,t 


if  (t+t )L  is  odd 

Tr(t2-t2) 


(C-7 ) 


(C-8) 


a 

S  -  b 


(C-9 ) 


Xr  =  ru  cos  (^)  n 


(C-10) 


it.  [(KM-k)2  +  X^l 

n  L  J 


(C-ll) 


v,  _  2kL  . 
br  -  ~T~  Xr 


(C-12) 


o  =  KM 

r>*  OT. 


2L  [Xr+tn/L]^  -  (KM)2  O^tt/l]2  -  (KM)‘ 


(C-13) 


r,t  2L 


[Xr+trr/L]2  -  (KM)2  [xr-t-/L]2  -  (KM)‘ 


(C-14) 


_  1 
2L 


Vtn/L 


\  -trr/L 


[xr+tn/L]2  -  (KM)2  [xr-tn/L]2  -  (KM)‘ 


(C-15) 
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<C-16) 


X+tn/L  X  -tn/L 

-1 _  4-  - - - 


r,t  2L  LCvWl]2  -  (KM)2  Cxr-tn/L]2  -  (KM)2J 

gr;t,t  =  er,ter,t  +  cr,tcr/t 

(C-17 ) 

hr,t,t  =  er,tcr,t  +  er,tcr,t 

(C-18) 

<1 

Et  =  21  ^ardr,t+brfr,t^ 
r=l 

(C-19) 

q. 

Et  =  21  ^arer,t+brcr,t^ 
r=l 

(C-20) 

<1  r 

Gt,l  -  Z  [(“rSr,t,trtrhr,t,t)  “s 
r=l 

+  ^hr,t,t+Vr,t,t>  sin  XrL] 

(C-21) 

q. 

®t,t  =  21  [argr,t,t  +  brhr,t,t] 
r=l 

( C -22 ) 

Procedure  1,  just  outlined,  is  suitable  only  for  small  values  of  m,  L,  and 
the  argument  of  the  Bessel  functions  Tu£  •  For  large  values  of  these  param¬ 
eters  the  number  of  terms  that  must  be  carried  in  the  series  approximations 
becomes  prohibitive.  Roughly  speaking  the  procedure  is  practical  for  m  ^  4  , 
L  £  2  ,  and  Tu§  £  10  . 


49 


Procedure  2 


From  the  definition  of  the  chordwise  deflection  function  given  in  Appendix  B, 
it  is  clear  that 

§mJ-Ja(l+0_C?)  =  Cm,J-Xrm(1+G-cP>  +  Dm,  > 
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where 


f  (v)  =  sin  y  v 
m  m 


sin  y** 
Tm 

sinh  v 
Tm 


sinh  Ymv 


It  is  therefore  convenient  to  introduce  the  following  general  integral 
definitions 


(9-p)f£(9)  dQ  = 


F^m(9-P,0,a) 


b 

a 


*  FS,m(b-p'b''l)  '  Fi 


m.m 


(a-p,a,a) 


where  from  [s] 


-  J_  -  4<a*3)el(n-p)49)<9> 

4vt  L m 

Ym 


L 

Tm  c 


Lpa*2)0](e.p)fi0)(6)  +  fiL(a+l)d(e.p)4«)(9) 


9  f[(<t*g)€(9.p)f(P9)(9) 
4  m  v  m 

+  (4c‘9>(9-e)438>(9)  + 

8y|  >• 


and 


(C-27) 


(C-28) 
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*fe,a(e-p,®,«)  -  -L_  .  ffats)e](e-p)fs(e)  - f<»e><9-p)fp><e) 

in  \  i 

+  f^a+1^(e-p)429^(0)  -  4(a+2)e](e-p)4e)(e) '  (c-29) 

w 

/  Q  ) 

In  these  definitions,  f^a  '(0)  for  a  =  0  is  fm(G )  ,  and  not  unity  as  the 
notation  seems  to  suggest.  Using  the  definitions  (C-27)  to  (C-29)  the  aero¬ 
dynamic  integral  becomes 

L  nl 

Im,m,u  =  -2-  /  Pu^"1+9)Gm  m  4(9)  dc? 

i-1  Jo 

L-l  r  1 

+  X  /  Pu(£-1+9)H_  (tp)  dcp  (C-30 ) 

1=1  Jo  ’ 


(l+e-qOfij^  j(e)  d9 


(l,<p,0)  -  FSjIn(l-<p,0,0)^ 


L 

"y~  c  -  c~  - 

-£— ,  m,j t-i  m,j l 

Z=Z+1 


m,m 


(-1, 1-9,0) 


L 

_Z  Cm,J-Z%J 

Z=Z+1 


F5,m(°^'°)  "  Fm 


m.m 


(-9,0,0) 


L 

Dm,  l-zFft,! 

£=jJ+l 


+  ^ 


Fm,m(9A,0)  -  FSjm(0, 1-9,0) 


(C-32) 

X=£+l 


integrals  f  PuU-l+<p)Gfi  _  ^(9)  d9  and  f  PuU-l+9)%>m^(<p)  d9  in 

;  *  Jn 


The  i 


(C-30)  are  evaluated  numerically. 
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APPENDIX  D 

EVALUATION  OF  INTEGRALS  Jm,m,Km,m,Rm,tn,  S,  T,  Q 


The  integrals  Jg  m,  Kg  m,  and  Rg  m  involve  the  chordwise  deflection 

functions  $  and  their  derivatives.  Therefore  introduce  local  coordinates 
m 

and  deflections  from  Appendix  B, 


0  =  -  -  (£-1)  (jC-l)a  sxila  (D-l) 

a 

~  ^  . L  (D-2) 


The  integrals  can  then  be  written  as 


n  L  L  nl 

/  5m(x)§fi(x)  dx  =  ae  (n-3) 

■  o  Jo-1  J  O 


Jm,m 


%,B  "  f  Z  [  4“,<0)W9)  ae)  <D-4) 

'  ^  n  i»  =  1  >-'n  ’ 


£=1  O 


%,»■/  4k)<x>5«<x)  dx  -  J  4!]<e>we)  de  (d-5> 


£=i  uo 


These  integrals  can  be  evaluated  conveniently  by  referring  to  the  Fg  m 
functions  introduced  in  Appendix  C,  Eqs.  (C-28)  and  (C-29).  There  results 


Jm,m  -  \  ’  FS,m<°'0'0)  [  51  < cm, iC«, «+D-,  A, t > 


1=1 


{FE,»<°'-1’0)  -  >a,m<1,0,0)]  Z  (»-6) 

L  J  ^=1 
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^m,m 


F-  (1,1,2)  -  P: 
m,mN  3  3  1 


m,m 


(0,0,2) 


L 

X  j&+r\n,£%,4) 

4=1 


+  I  Fm,m^0,"1,2)  -Fm,m(l'°'2) 


Z  I'iAA.Ai1  (D_7) 

4=1 


Sm,m  =  [Fffi,m(l,1,l)  "  Fm,m(°,0,l)  ]  ^  (Cm,£Cm,£-Dm,  4Dm,4) 


+ 


(o,-i,i) 


-  F_  (1,0,1) 
m,mv  *  *  7 


L 

( cm,  4Dm,  jj-Dm,  jJCm,  4  ) 
4=1 


(D-8) 


The  integrals  S,  T,  and  Q  involve  the  spanwise  deflection  shape 
and  its  derivatives.  Introduce  local  coordinates  and  deflection 


T]  =  y  -  (k-1)  k-Hysk 


(D-9) 


*n  =  *n,k 


k  =  1,2,3, 


(D-10) 


and  subsequently  obtain 


pN  „  N  pi 

-  J  *>-Z  / 

k=l  Jo  ' 


(D-ll) 


T-  I  iZy)\  ay  -  Z  I  «1 


n  ri 


k=l 


t(2D)/ 


(D-12) 


J  o  k=l  “  o 


(D-13) 
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Since  these  integrals  depend  on  the  spanwise  deflection  shape ,  it  is  again 
necessary  to  consider  each  of  the  five  physical  situations  of  interest  to  the 
present  analysis.  These  situations  are  listed  in  detail  in  the  report  body 
and  in  Appendix  E.  It  turns  out  that  for  cases  1,  2,  and  3 


N 

s  =  ^  — 


k=l 


*  [b„,ic431,>'i>  -  In,k45T|)(<»] 

* 

-  -  Vk^’H  [5n,k42’’)<1>  * 

*  [ffn,kfi,|,(0>  -  5n,k?n1,,<1)]  [^k^0)  *  Vk5^1)]  ]  <D-16> 


N 


T  -  Z  -ij  -v‘[5n,k4T|,<1>  -  5„,kfn’1,(°)]  -  -  0n,k431,)((>)]‘ 

k=l  4v  L 
Tn 


Kk^’t1)  +  \k*ia|)«»]  [®n(kIiS1|)t1)  -»n,k431,)(°)] 

[En,kir)«»  +  Sn,kIn2ll>(1,][ei>,k4Sn>(0)  ’ 


(D-17 ) 


-4 

Q  =  v  s 

ln 

For  cases  4  and  5  these  quantities  reduce  to  simply 


(D-18) 


S  =  1 

T  =  Q  =  0  (D-19) 
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APPENDIX  E 

EVALUATION  OF  FOURIER  COEFFICIENTS  Bn;U  AND  THE  INTEGRAL  F(u) 


The  spanwise  deflection  shape  ijf  (y)  is  expanded  in  a  sine  series 


E 

u 


n.u 


sin 


UTT0 

B 


(E-l) 


where  <f)  is  the  transformed  coordinate  (Pig.  3) 


<f>  =  y  +  x 

B  =  2\  +  N  (E-2) 


and  x  is  a  constant  to  "be  evaluated.  The  coefficients  Bn  u  are  given  by 


B 


n,u 


£ 

B 


(fn(0)  sin 


d0 

B 


(E-3) 


These  coefficients  are  to  be  evaluated  for  each  of  the  five  physical  situations 
to  be  analyzed.  The  situations  are  repeated  here  for  completeness. 

(1)  A  panel  array  as  shown  in  Fig.  1  with  a  finite  number  of  panels 
in  both  the  chordwise  and  spanwise  direction.  The  array  is  bordered  by  an 
inflexible  surface  extending  to  infinity  in  all  directions. 

(2)  The  same  as  (l)  except  that  the  array  extends  to  infinity  in 
the  spanwise  direction. 

(3)  The  same  as  (l)  except  that  the  array  is  flanked  on  the  sides 
by  vertical  walls  (a  wind  tunnel  installation  for  instance)  an  arbitrary  dis¬ 
tance  from  the  edge  of  the  array. 

(4)  The  same  as  (l)  except  that  the  array  has  one  panel  in  the  span- 
wise  direction  (N=l)  and  the  side  edges  are  free. 

(5)  The  same  as  (4)  except  that  the  array  is  flanked  by  vertical 
walls  as  described  in  (3). 
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Before  evaluating  the  Bn  u  it  is  useful  to  first  determine  X  for  each 
case.  Geometrical  aspects  of  the  expansion  (E-l)  are  shown  in  Fig.  3  for 
cases  (1)  and  (4).  The  finite  span  deflection  shape  is  depicted  in  Fig.  3a, 
and  the  periodic  expansion  of  shown  in  Fig.  3b.  The  domain  of  aero¬ 

dynamic  influence  is  illustrated  in  Fig.  3c.  All  points  inside  the  shaded 
region  are  aero dynamic ally  coupled  to  the  panel  configuration  being  analyzed. 
It  is  seen  that  for  X  £  sL/2M  the  periodic  reflections  of  are  outside 

the  coupled  region  and  the  effect  of  an  isolated  panel  array  is  obtained. 

For  cases  (3)  and  (5)  the  effect  of  vertical  walls  is  obtained  by 
introducing  fictitious  image  panels  on  each  side  of  the  real  panel.  The  geo¬ 
metrical  aspects  of  this  situation  are  shown  in  Fig.  4.  The  wall  locations 
must  be  planes  of  symmetry  in  the  flow  and  this  situation  is  realized  when 


—  N+2A 
M 


The  desired  aerodynamic  isolation  is  obtained  by  taking  X  =  N+3A/2  . 

The  case  of  a  panel  array  extending  to  infinity  in  the  spanwise' 
direction  (case  2)  offers  some  difficulty.  In  accordance  with  the  discussion 
in  the  report  body,  is  restricted  to  cases  in  which  the  deflection  in 
adjacent  panels  is  antisymmetric  in  which  case  X  =  0  and  B  =  1  . 

In  summary,  then,  X  is  taken  as  follows. 


Case 


1 

2 

3 

4 

5 


X 

sL/2M 

0 

larger  of  sL/2M  and  N+3A/2 
sL/2M 

larger  of  sL/2M  and  N+3A/2 


The  coefficients  Bn  u  are  evaluated  in  the  following  sections. 
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Case  1  -  Panel  Array  of  Finite  Spanwise  Extent 
The  deflection  is  given  by 


0 

0  £  0  £  X 

X  £  0  £  X+N 

ft  =  0 
vn 

X+N  £  0  <:  2x+N 

(E-4) 

From  (E-3)  there  follows 

n  X+N 
■  1 1 

(0)  sin  d0 

B 

(E-5) 

Introducing  local  coordinates  and  a  local  deflection 

I]  =  y-(k-l)  ; 

k-1  £  y  £  k 

k  *  1,2  ....  N 

(E-6) 

yields 
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B. 


3  N 

2B  £ 


'  t.4-4  „4  4  v,  | 

B  v  -u  tt  k=l  L 
Tn 


ea,k?£sn><i>  -  \,kf»311>(o) 


-(r)2[Sn,kfn1')<1>  -  “n,*^0)] 


sin  (k+x) 
B 


•  (?)  Kk^V)  +  Dn,kfi21l)(0)]  cos  f  (k+X)  -  {  Cn>k‘43Tl)(0) 


-  Dn,kfi31l)(l)  ~(A^\o)  -5a,k41l)(l) 


(?)[Cn,k427l)(0)  +  ^n,k421l)(1)]  cos  ?  (k+X-l) 


sin  —  (k+X-l) 
B 


(E-7) 


where  Cn  k  ,  Dn  k  and  fn(T])  are  defined  in  Appendix  B 


Case  2  -  Panel  Array  of  Infinite  Spanwise  Extent 


The  coefficients  for  this  case  are  obtained  from  (E-7)  by  setting 
k  =  N  =  B  =  1.  There  follows 


2uit 


n,u  _4  .  .4 

Yn-(^TT) 


(E-8) 
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Case  3  -  Panel  Array  of  Finite  Spanwise  Extent  with  Vertical  Side  Walls 


The  deflection  is  given  by  (Fig.  4) 


*n=  0 

0  <■  \ i  <■  X.-N-A 

V  = 

X-N-A  S  0  <:  X-A 

♦n“  0 

X-A  ^  0  ^  X 

X  £  0  £  X+N 

=  0 

X+N  ^  0  s  X"*"N+A 

=  «n 

X+N+A  S  0  S  X+2N+A 

*n  =  ° 

X+2N+A  £  0  £  2X+N 

From  (E-3)  there  follows 


+ 


2 

B 


\H-2N-fA 


U0)  sin  m.  ^ 

11  T> 


which  gives 


(E-9) 


(E-10) 
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B 


2B" 


N 

'n,u  44  4  4 

B  V  "u  TT  ^-1 
Tn 


C  „f(3Tl)<l)  -*  J<31,)(0) 

n,k  n  v  7  n,k  n  v  / 


-  /Hlf  fc  .  f(71,)(l)  -  D„  ^(O)!  [sin  HU  (k+X-N-A)  +  sin  HE  (k+\) 
\B  )  L  n>k  n,k  n  J  J  L  B  B 


+  sin  HU  (k+x+A+W) 
B 


"  [cn,kfi2ll)(1)  +  Dn,^211^0)]  [cos  (k+X-N-A)  +  cos  HE  (k+x) 


+  cos  HEL  (k+x+A+N) 
B 


1  _(3T|)  _  .(311),  x  utt 

C  .  f  (0)  -  D  ,  f  (1)  -  — 
n,k  n  v  7  n,k  n  v  7  3 


C  -f^(0)  -  D  f^(l) 
L  n,k  n  v  7  n,k  n  v  ' 


X  [sin  HEL  (k-l+\-N-A)  +  sin  HU  (k-l+x)  +  sin  HEL  (k-l+X+A+N) 

Lb  b  b  j 


UTT 


K^C)  * 


x  [cos  HE  (k-l+X-N-A)  +  cos  HE  (k-l+X)  +  cos  HE  (k-l+X+A+N)  f  (E-ll) 

Lb  b  b  j 


Case  4  -  Panel  Array  with  N  =  1  and  Side  Edges  Free 
The  deflection  is  given  by 


♦n  -  o  0  £  0  s  X 

♦n  '  1  *  5  *  f 

1fn  =0  X+l  *  0  £  2X+1  (E-12) 
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From  (E-3)  there  follows 


Bn  =  i  f  sin  HE  (0+X)  d0  (E-13) 

n>u  B  J  Q  B 


or 

B  =  ?_  fl-(-l)Ul  cos  HE  X  (E-14) 

n>U  UTT  *■  ■*  B 


Case  5  -  Panel  Array  with  N  =  1  ,  Side  Edges  Free  and  with  Vertical  Side  Walls 

This  case  is  a  special  form  of  case  3  with  k  =  N  =  =  1  . 

Straightforward  integration  gives 

B  =  —  I  cos  HE  (x-A-l)  -  cos  HE  (x-A)  +  cos  HE  (x) 

n'u  UTT  [  B  B  B 


-  cos  HE  (x+l)  +  cos  HE  (x-i-A+l) 
B  B 


-  cos  HE  (X+A+2)  l  (E-15) 

B  J 


A  second  purpose  of  this  Appendix  is  to  evaluate  the  integral 


f  rk 

F(u)  =  /  sin  HUE  d0  (E-16) 

J\  B 

introduced  in  Eq.  (29)  of  the  section  dealing  with  the  aerodynamic  terms.  In 
evaluating  F(u)  it  is  again  necessary  to  be  cognizant  of  the  five  physical 
cases  of  interest.  Comparison  of  (E-16)  with  (E-3)  shows  that  for  cases  (l), 
(2),  and  (4) 

F(u)  =  BBnjU/2  (E-17) 
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For  cases  (3)  and.  (5),  F(u)  can 


be  evaluated  in  closed  form  to  give: 


Case  (5) 

F(u) 


Case  (5) 


4  N 

_ if - £ 

(By  )4"(u1t)4  k=1 


r  f(31l)d)  -  d  -f(3^(0) 

cn,lrn  un,k  n  '  ' 


sin  —  (k+\) 
B 


-C  vf(31>>(0) 

n,k  n 


+  D  ^>(1) 
n,k  n 


sin  HU  (k-l+X) 
B 


HH 

B 


HU 

B 


*  Sn,.42T1>« 


cos  —  (k+x) 
B 


cos  —  (k-l+X ) 
B 


(E-18) 


(E-19) 
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"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof:’ 

— National  Aeronautics  and  Space  Act  or  1958 


NASA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 


TECHNICAL  REPORTS:  Scientific  and  technical  information  considered 
important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

TECHNICAL  NOTES:  Information  less  broad  in  scope  but  nevertheless 
of  importance  as  a  contribution  to  existing  knowledge. 

TECHNICAL  MEMORANDUMS:  Information  receiving  limited  distri¬ 
bution  because  of  preliminary  data,  security  classification,  or  other  reasons, 

CONTRACTOR  REPORTS:  Technical  information  generated  in  con¬ 
nection  with  a  NASA  contract  or  grant  and  released  under  NASA  auspices. 

TECHNICAL  TRANSLATIONS:  Information  published  in  a  foreign 
language  considered  to  merit  NASA  distribution  in  English. 

TECHNICAL  REPRINTS:  Information  derived  from  NASA  activities 
and  initially  published  in  the  form  of  journal  articles. 

SPECIAL  PUBLICATIONS:  Information  derived  from  or  of  value  to 
NASA  activities  but  not  necessarily  reporting  the  results  of  individual 
NASA-programmed  scientific  efforts.  Publications  include  conference 
proceedings,  monographs,  data  compilations,  handbooks,  sourcebooks, 
and  special  bibliographies. 


Details  on  the  availability  of  these  publications  may  be  obtained  from: 


SCIENTIFIC  AND  TECHNICAL  INFORMATION  DIVISION 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

Washinston,  D.C.  20546 


